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:  R  I 

INTRODUCTI 

Solid      -    jr  powder  blending  is  an  operation  by  which 
solid  materials  in  particulate  forr  are  scattered  rar.      ?r   unifor 
aniong  each  other  by  the  movement  of  the  particles  in  a  mixer.   '.'  -      ..- 
tion  may  be  broadly  defined  as  anv  process  that  tends  to  elir.ina 
inhomogeneities ,  or  to  reduce  existing  gradients.   On  the  other 

ty  be  performed  to  achieve  an  acceptabl- 
control  rates  ol     t  transfer,  mass  transfer  or  chemical  reaction. 

ziough  it  is  one  of  tne  oldest  unit  operations,  it  is   til]  one  of  the 
most  wid._     :>loyed  in  industries.   For  example,  it  is  essential  in 
plastic  processing,  ore  smelting,  fer-      r  production,  pharmaceutical 
Lon,  food  manuf.      and  catalytical  t  is  of  chemical 

in  spite  of  the  practical  importance  of  I 
operation,  the  fundamental  theories  underlying    .s  ope-:      and  the 
procedure  for  apprai  .         .  formancc-         s  have  not  I         -y 
est.: 

Solids  mixing  pr^  •  indu-      are  mult, 

in  •  -ol ended,  and  t:   . 

of  utmost  practical         mce. 
Ch  literature  on  mixing  and  m .      -.  of  raulticom- 

.  r  t  i  ]  . 


The  concept  of  the  contact  number  is  extended  to  a  inulti component  solids 
re  In  the  completely  mixed  state  In  chapter  HI.   Such  an  approach  has 
been  applieo      to  the  analysis  of  a  binary  mixture  (1].   Chap:      |  also 

an  estimation  of  the  mean  cont..  spot  sa         j    the 

re  -tulation  of  mixtures  at  various  concentrations  of 

the  k>.      ponent.   The  number  of  spot  samples  necessary  for  a  specified 
value  of  the  relative  standard  error  is  derived  in  this  chapter. 

In  Chapter  IV,  a  Dirichlc      .nomial  model  is  proposed  to  describe 
y.ture  in  an  incompletely  mixed  stau  .       ates  c   I 
parameters  of  the  Dirichlet  distribution  are  derived  and  the  precision  of 
the  estimated  population  mean  contact  number  is  assessed.   This  mod      <nds 
the  beta-binomial  model  for  ah;        ture  in  an  incomplett 

state      to  the  multicomponent  case.   Computer  simulations  are  included 
to  demonstrate  the  effectiveness  cf  the  derived  model. 

To  characterize  the  homogeneity  of  a  solids  fixture,  a  mixing  index 
usual 1      .  oyed  as  a  measure  of  the  degree  of  mixedness.   Over  th: 
rent  c       a  for  the  degree  of  mixedness  have  bet.      I  ewed  and 

Fan  et  al.  (  3  ].        unately,  the  general  criterion  for 
^ompor.  i?  still  lacking.   Chap*. 

presents  some  approaches  to  define  the  mixing  indices  for  rail t icomponent 
. 
ipter  cant  conclusions  and  proposes  several 

pror      areas  of  future  research. 
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Numerous    pul  Lications    or.   binary    sclics    mixtures    ar.c  Lnar; 

systetos  are  contained  ~r.  a  number  of  surveys  and  review  am:les.  including 
those  published  by  Ueidenbaum  [1],  Bourne  _-],  Far.  et  al .  [3-6],  Yalr:  tin 
[7-9],  Yar.:  [10],  Williaa  [11,12],  Cooke  et  al  .  [13].  Since  the  Study  on 
mixing  of  binary  panicle  systems  and  their  oharacterioa  t  loo  is  essential 
to  the  understanding  of  multicomponent  systems,  the  recent  publications  on 
binary  mixtures  arid  nixing,  vhicb  are  not  included  in  the  previous  survey 
and   review   articles,    are   alsc    reviewed    in    this    article. 

There  -are    two    broad    and    equally    important    aspects    in    the    stud;  Ldl 

.r.g.      One    is    concerned   with    the   characterization  of   a   mixture    and    the 
other   wit  nechanisit  and   rate  by  which    the    state    of    a   mixture    changes • 

Loth    aspects    are    included    in    this    review. 

1.  m:e:;t  PUBLICA  v  solids  MIXTURES  ANI  MIXIH 

7 section  rief  review  I   the  reel  res 

sclics  particle  systems. 
1.1.        :es 

A:  '-eld  [1^],  three  crit   .   are  I 

Tnese  are     -  cn-sar.;.  I  Intensity  of 

segreg  .  on  (s:        .  c  str- 

and Kersey  [15]  utilized  Suslik's  concept  of  homog 
to  define  the  homogeneity  of  a  \  ,  which  has  zero  standard 

ion  for  the  cone      :on  of  a  cr.inor  ingredient  be:-     samples. 
.  senhar t-Rothe  and  Peschl  [16     Eroduced  concepts  of  testing  equipme 
for  bulk  materials  har.      and  also  described  tests  for  characterising  the 
flow  propertiei      >wders .   tfickl  [17]  blended  I 
particle  size  t  utions  to  achieve  a  target  particle  size  distribution. 
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|    developed   a   new  for  det  ;le   size 

-    < 

are  two  approaches  to  the  determination  of  the  nixing  index  based 

on  the  contact  number.   One  is  based      ■•  coor      .jn  number  sampling, 

Chi  sp         ng.   By  these  two  approaches,  Shindo  et  al. 

posed  th<.         models  of  the         ition  cf  the  conta 

an  incomple:         I  state. 

. 

Ra:  beds.   The 

ef:  :atior.  or.  the  hydrodynamics  of  gi  |   Flow 

necessary  : gence  of  a  reg      :  article  stream 

wert      :ssed.   He  alsc  showed  that  the  flow  pattern  in  the  .        the 

EunctiOl  parameters. 

.1]  re.         ree  mechanisms  cf  segregatic  r . :   :: 

segre.  I  particles,  and  the  rise  cf    coarse  particles 

KoseV.i  e:  al.         tdiod  segregating  phenomena  in  a  ,-.:::- 

.  arum  mixer.   Steedman  e:  il.  [23         -gated  sclids  segre- 

a  coal  gasification  burner.   Sugimoto  [24]  Itudlod  tolidl  Mixing 

frc      ■  segregation  characteristic:      .xrures. 

urega:         -  lex  par:  tier..  a 

phenomenon  caused  by  a  shifting  of  fines  toward  the  bottom  of  the  container, 

has  been  reported  for  the  systems  with  particle  diame:        Cht  range  of 

:>].   Parson  (26]  extended  thJ 

particle  I        -.    and  without  the  presence  of  agglomerates. 

Row  a:  revieve:  ~k.s 

on  z-z    se,-  gas    :'  red   beds.      7.\ey    sis:    reported 

lOOt    imi    I  eroble: 
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Burgess  et  al.  [28]  usee  zas    : laicized  bees  to  study  solids  mixing 

ar.d  segregation.   They  found  that  mixing  commences  a:  gas  velocities  just 

above  the  minimus  fluidization  velocity  of  the  non-segregat ir.g  component 
for  particle  systems  of  equal  density,  but  nixing  does  ret  commence  until 
the  gas  velocity  exceeds  the  minimus  fluidization  velocitj  of  the  se. 

component  for  systems  of  different  density.   They  aisc  developed  a  model  for 
segregation  and  r.ixing,  which  is  capable  of  predicting  the  cc      behavior 
found  experimentally . 

-tzgeraid  et  al.  [29]    used  a  fiuidized  bed    to  stuc;.  solid       Dg 
and  calculated  solids  flov  patterns  by  use  of  a  general  ;      Eaentsl  moot. 
Pvengarajar.  et  al .  [30]  employed  the  single  phase  backflow  cell  mo^ 
simulate  the  solids  mixing  ir.  a  fiuidized  bed  coal  combustor. 

-y  and  Merrick  [21]  developed  a  mathematical  model  to  C     - 
eral  concentration  gradients  of  reacrant  solids  ir.  a  1  | 
to  predict  the  effect  of  coal  feed  spacing  on  the  com'r  .   a 

bed   boiler.   Lateral-solids  mixing  diffusivltiei 

in  a  5  ft.  diameter  bed. 

economic  advantages  of  using  a  continuous  process  .        - 
accepted.   However,  processes  involving  the  h.  »C easing 

particulate  solids  are  often  designed  as  ba:      ocesses.   Such  systems 
should  be  examined  to  determine       would  be  advantageous  to  I        9n- 
tinuous  nixing.   Williams  [11  r.uous  mixing  of  solids, 

in  which  the  advantages  of  continu  <  equipment  for  continuous 

•  are  also  report!     :yanami  et  al .  [32]  developed  a  continuous  i 
process  for  fine  and  cohesive  powders.   Th<      .cal  cylind'     it  provided 
with  an  impeller  of  multistage  comblike  blades  has  bet-     .en  effective. 


- 

Carstensi  -  ed    that    the  otxing  of  a  b:  ;>n- 

•pherical   particles    -  rough  surfaces   obeys   tc  rial   equa:  as 

long  as    the  me.  the   :-  actions  are    id-  and 

Hwang  Lng   of   a  racer 

Iht  1  exper  i- 

■i  analyzed    .  ;ns  of   a  diffusion  ao  Lt   if    assumed 

■a    function  of    the  a: 

stream. 

e  t  e  d  I  salk 

doc  f  grains  in  a  notior.Ies;  .rer    (K. 

»r.      ~  pcsec   a   discrete   steac; 

-     •  segre;  is 

mod    .  degree  ecnes' 

the  -ites    (  ^e   particle-  blended   by  :he 

par •  tr. 

Investigated  the  Mixing  proces 
bat  -  |  35  ]    usee   a    p: 

1  •    .  Bg  and    her  .jir.g  ;ls. 

-ater  h  ed  funds  Lssm .     Yaaagucl 

[AO]  developed  so=e  rate  equations  for  sclics      g. 
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3.      MDLTK  EN1  MIXTl  CNG 

rne  publications   or.  Bulticom]  >i  et  c  nixtures  a-  review 

ansively    in    this    section. 
2..1  res 

Gayle   e:   ai.     »-l  ,   used    chi-square    as    a    cr.ter.  the 

.:'    a   mult  icomponent    solid]     mixture    inst  of    a    c 

-r.g    index  based   or.    the    sample   variance    (or    sta:  :ion). 

Ace;  •  tc    ther., 

,    n  (0.  -  E   • 

■i     *      i  o 

i-1    Ei 

whe 

■  mi  --s, 

.  -   uare  v  :  :  freedom, 

n  & 

=  observed  number  of  pa  Lee   of  any  given  color  in  the  L- 
sampie  , 

■  corresponding  I  in  the  i-th  sat.;.   .  .  asec  oi 
average  distribution  of  components  in  the  mixture. 

They  also  proposed  a  segregation  index  as 


-  -? — h 

s    2    2 

- 


I   -  numerical  index  which  indicates  degree  of  segregate      the 
>.  •  obaer       -  *qu«r«  re, 
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-square  :     i  randoc      re, 

-  e>.:  i-»qu-         :e  segregate        •• 

Be  loped  a  formul..         variance  of  a 

.ven  p.  samples  from  a  granular  material 

a  k-  -  .on  as 


_ 
q  (1  -           v  -  v  ) 
.  -J — t = 5 (3) 


vht 

icb  samr.  . 
:T.e  per  V.  btJ 

III  aggregate  w'r. 

random  samples 

t  volume  V  and 
«  average  volume  cf  a  particle  in 

as 
v-      -  ,2v2  -  ...  -  q. 


whe:         q  m  fractions  of  the  part 

1 


aggregate . 
Stange  [43]  characterized  the  quality  of  «".      :re  consisting  of 

s  of        cnt  sizes,  denoted  b    ...},...  He  developed  the 
Cbl  variance  of  one  of  the  component  g  the 

component?  as  z.  us  compo: 


:i  -  b 


r  i  -  Pi 


V — ~   7.   Ef  «  ),+  P,  (If  w  ),+  P.   I:  •-.) 

1    :  :  -   2    j  j  2 


-  etc 


«0 


wnere 

7.,  p  etc.  =  «  Lent  proportions  of  cocpor.     .  2, etc  . 

(If.w.J.     =  effective  mean  particle  weight  of  co-.oonent  i. 

expression  ca  extending  Bus  111     .:..  (3)  I    over 

Dulti component  system  .---■ 

Knc::  .-:.  ir.i^lycec  £  mixture  of  particles  of  ;i:es 

mg  a  technique  used  Li  the  study  .  -  ,  : 

various  types  in  a  renewal  process  up  t   a       time.   Be  del 
for  the  i»yi  |   Lc  means,  variances,  and  covariances  of  t         I  Ol 
rent  types  of  particles  as  the  total  \ 
.  ed  asymptot.  .  is 


r 

.   Vh  v   (i)  L 

?-  «   i 

3 


. 


.'--.':■-. 


-  u2(i)J  ,        . 
■ 7, =■  [w  -  2r.,-(i)j 


where 


h   •  number  fraction  of  p.:     ••  of  I 
\   ■  volume  of  particles  li     h  sample, 

an  of  the  ■!*<  on  of  particles  «. 

m(i)  -  second  moment  about  the  Loo  of 

particles  c: 
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t    ■•   .  .  -■  s  . 

Sc  .plied  the  :  th«  as      ic 

to 
:.s ,  and  to  c  ..ge 

..nee  and      .ance  are  et 
|  roposed  I  'a 

assu:-        :he  popul   .  covariance 

are     -  vai  pi  pesed  as 


I  s  , 

is  asympotc I .  • ed  as  a  aultivariai 

,  1  •  .Dutior       -  »  asympoi 

uni  oroug:     Lxad  aggregate.   1      .d  be  close  to 

tnt  Sggragai  .  -nd 

•  cd  or.  Statistic  does  not 

a  gooc 
a  very  1. 

Bu  duced  a  co:  jgene-' 

.  .  .  .  ,  a 

s : .  .       .  •  .   : 
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is  difir.ec"  as  the  negative  logarithn  :f  the  sample  veight. 

•  .50.  exter     this  concept  to  e  nulticomponent  mixture  by  i 

the  standard  deviation  of  an  ingredient  1  in  a  cor        ra 

fixture,  eqn.  (^).  They  assumed  that  the  components  must  be  re. 

san>u  par tii-tr- size  level   prior  tC  mixing,  i.e.,  particles  of  eat.     .one: 

have  the  sa:.-  particle  weight. 


(f.vO.  =  E(f> 

3   J  -     :  : 


=  V 


re  H  is  the  effective  near,  weight  of  all  particles  in  tne  fixture.   Tnus , 

eqn.  (i)  reduces  to 


::  =  vV;-v 


2 

:.    C\  ■    0.01    and   U    =   W.    gives    a    .  rent    equaticr.    :;r    thl    homogeneity    of    ar 

redient    in   t  .Disponent   nixture   as 


j)    -   -log   [p. (100   -    p   )    wj     ,      j   •  1,2,. 


(9) 


is  can  be  used  fa  single 

component  in  a  mixing  operation  or  to  compare  thl  of 

powders  necessary  f      *ng  th<  horaogent      .a 

.  componer.        operation. 

Cook  and  evaluated  a  Kauta  DX  (  g  of 

a  multicomponent  tablet  preblend.   Four  ingredients  (phenobarbitone, 
butobarbitone,  quinalbabitone  and  lactose)  were  i  :ng 

operation  samples  were  removed  fro  and  ass.<      r  each  ingred: 
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gas    1  -ornate.  11  ••   of   i  l nt 

as    cal  ed  accorc  (9).       Eac 

beh  .e  Banner   during   siring.      The   control  jae   is 

nee  a:   ary   one   tioe  one 
segre.  oes   not  meet   the 

jmoge: 
[53]    proposed   a   discrt  model 

of   a   i  omponent    homogeneous   p. 

*  r  r.      The    cor.  .ance,    -    ,    of    the  .it 


;■  -  <vJ 


-    -re 


C.  (K)  -  number  concentration  cf  component  J  ::.   • 
(C/'b  ■  ec  conct 

■  n  u- :  t  r  f  ractioi 

number  cf  ceils  in  a  mixer . 

r  components  in  a  mixer. 

le  composition     tht  completely  segregated  state, 

.    was  •    c   model   as 

o 

*         GvJ        L'-       t0>"1]2 
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where 

P       =   transition  probability   of   the  particles    fr<  LI  i  tc    ze'. 

i£ 


C     (0) 

Q  .  .(0)  =  :.  ,  fj- 


The  variar.ee  reduction  ratio,'-  >  car.  be  vritten  as 

I  -  ^5  (12) 

o 

The  degree  of  rtixedness  r.ay  be  defined  as 

M  -  1  -  1  (13) 

La_  and  Far.  [53]  proposed  the  idea  cf  using  the  e       c:    a  sclic  ~- 
as  a  criterion  for  characterizing  the  mixture.  Ihe  entropy  of  the  sc  -  . 

r.ixture  a:  cell  i  car.  be  defined  as 

k 

5  (N)  ■  -   I   C.  .(N)£n  C   (N)  (1- 

1         j.l    J         J 

and  •    total  entropy  of  the  syster.  is  then 

a 
S(N)  -  ::..£.  (N;  L5) 

1*1 

the  systen  is  ir.  the  segregated         t  the  beg: 

process, 

£,(0)  -  0  .6) 

and  therefore, 

a 

S(0)  -   :   S .  -  0  (17> 

i 

i-1 
Th<    •  .  jdv  of  a  1  cha:.        .:ero  in  the  conoletelv  secre?ated 

state  to  a  constant  val        ■  comole:      .»ndoti>  state.   Hen..  .  Che  en::. 
of  a  i  can  be  sin.il;.-       ned  as  a  measure  of  the  decree  I        ess 

of  a  i     c. 
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I    a 
mean  vectors   fti 
chat  .nee  a.*  «    st*:  I 

-  -ies   has 
bee:  ited.  :-e   degree 

e  at 

:     -    . 

*  --^ S) 

-r 

re 

-  ... 

■  thf    cow.  segregated 

:e, 

I    |   »  <  t  of    th«  the  co-  •    • 

- 

ihe 

■pi  1 1  (  •.-.-.-• 

tveral   n.cnosizec  .s,    Soroer    [55]    a< 

1  usm  .    •     , 

Lzt  -    t . , 

3   -    .     .     .     .  ....... 

-  :.ie   greatest    size   anc  e 

I    assumption,    the    f c I  can  be   used    to    calculate 

.ncc   of    the    ■  I  Lc    horoogt  as 

j-l 
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c   =  — , 

: 

1  • « q* 

X.  =  volume  concentration  in  a  sample, 
q  =  theoretic  volume  concentration  of  component  i. 
3.2.   Hixing 

Gayle  et  al .        arried  out  nixing  experiments  b-        1  mixil 
Ich  was  designed  for  mixing  samples  of  coal  and  coke.   Gra-  . . 
tch  were  .      r  in  surface  characteristics,  shape  and  densi:   I  at  iif- 
I  .  .  vi  r«  used.   During  nixing,  the  value  of  ch i-scuare,  as 
.;.  ecr..  (2),  decreases  with  an  increase  ir.  the  number  of  revolu- 
tion!        ---:ir.g  wheel  and  approaches  a  mean  lower  limit  asymptotical] 

2 
Tbf  initial  value,  x      »  -or  chi-square  (before  mi:  Lr    is  equal  : 

number  of  degrees  of  freedom  multiplied  by  the  number  of  par:.  Lei 

foi  each  temple.   The  number  of  degrees  of  freedom  is  defined  as  : 

of  .     of  data.   The  segregation  index  of  eqn.  (3),  v 

unity  to  zero  as  mixing  proceeds  from  co_;le:e  segregat.cn  to  cor. 

randomization,  was  used  to  establish  a  conventional  rate  equation  as 

dI^      e 
n  i 


where- 


dr       s 

I   ■  di-      if  segregation 

■  num.-  ons  o:        er, 

n  ■  rate  consta: 
8   •  consta: 
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La.     "an  (53]  developed  a  stochastic  node!             .,  of  ic-:.- 

cor        .mogeneous  :  a  s>e  |          >.er 

inentaliy  dc  ^abilities  of  a  I 

homogeneous  pa:          km.   The  oodel  pe:                  ihe  concen- 

and  the  degree  c:  .  copone     -„o- 
ger.                   e. 

et  al.  [54]  employed  a  component  solids 

•  s  hav.     LtntlCl             except  color 

I  used.  .  ance  mac: 

reases  as  rease? . 

is  define  -qn.  (18'.    reases  and  approaches  1  as 

the              ;  reases.  I   a  Milt     poner.  t  heterogeneous  z.. 

expe:              reases  up  and  decreas- 

asyr  to  tidily  to  1  tate  because  of  the  segre:       tenet- - 


-.     : 

review  of    th(  re   or.   binary   fixtures    a-  I  es 

r.ary    c or. poner. ts    systems  .  ized   bee  exter. 

the    last    several    years.      Zhll    tfl    a    net  onsequence   cf    the 

zed   beds    in  cany    areas   cf  material   process 
Lag   rar .  This    rev:ev    alsc    indicates    that    segrep  OX 

demixing,    which    is    the   counterpart   of  mixing,  than  nas 

recor.e    the   esphasis   of    the   re  .gatic 

sampling  of  toulti  M    and 

.ons    have  been   of    intense    interest.      Vet,    the    literature 

■ ance,    the    c:  ting 
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the  homogeneity  of  a  heterogeneous  Bulticomponent  nlxture  is  still  la 
Sir.ce  any  knowledge  of  Bulticomponent  sclids  mixing  and  mixtures  shouJ 
reduce  to  that  of  binary  mi:  ng,  tc  establish  a  bread,  systematic  ar.c 

analytic  approach  tc  this  subject  is  of  importance  and  significance. 
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obser        ar  of  particles  of  n:    .  to  color  in  t  I  L-l 


.on  probability  of  the  particles  from  eel. 
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q  ooretical  volune  concentration  of  component  i 
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1        qi 

Q,.       a.  ^ 

^  e  r 

.  nee  m 


II  -  18 

7        testing  statistic  as  -  ..  (6) 

t         required  tir.e  for  one  pass 

volume  of  particles  ir.  each  sar:.; 

vcl  mne  c  :"  a  p  - 1  "  i 

v         average  volu      a  particle  ir.  the  -ixture 

weight  of  a  sample 
W         weight  of  sample  necessary  tc  give  a  standard  de. -  ... 

ear.  weight  of  all  particles  in  : 

volur.e  concentration  ir.  a  sar.'.cla 

Greek 

u        total  nuDoer  of  cells  in  a  r.ixer 

rate  constant 

X  ■    anour.t  of  a  sar.ple 

•-(i;  lor  volune)  distribution  of  1 

as  defined  in  eqn.  (5)  = 

i-1  X 
II.  n<.  fraction  of  par-.      in  cell  i 

-covariance  natrix  of  a  xed 

:te 

variance-c>  trix   of   a 

o  indard  n   of   v  .on  of  Jom  samples 

of    constant   vo . 
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I  is  tak-.  :e: 

of  all  particles  which  are  in  contact  -.     :s  particular  particle 

er  der.    I       :*,  ar)d    the 
:alled  the  I  Such         5  is  called 

on  number  M     |  of  I      *.   For  exar 
D  number  of  the  t      ensior.     -  ic  packing 
•it  of  the  hexagonal  packing  I      -Tient  is  6. 
Let  A.  particles  be  key  particles  in  a  E      omponent  mixture 


con 


cles  :      same  size,  Ar  ,.-. I    .  The  number 

u  1      - . 


.es  of  conponent  A.  in  contact  vitta  a  key  p«  .3  def I 

as  the  contact  \ ributed  by  conponent  k.%    and  is  denoted 

1(0)'            e'  We  S6e  :"  r:C          1(0)  "  1'  C'>(0)  "  3' 

3.              of  the  Contact  Number 

I  complet  .xed  state,  particles  of  ea:r 

•^cted  to  be  randor  ibuted  in  a  sample  f]      e  cocr     .on 

number  sampling  of  sice  n*.           ,  the  distribution  cf  the 

contact  numbers,  C.  ,-.  ,  C.,, .,...,  O^ix/os  •                      •- 
Dution,  i.e. 


*0   L0(0  •     Cl(0) 
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Fig.  1.   Illustration  of  a  coordination  number  sampl . 
size  d*"( 
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conce  ~onent 

,C/rT  -.us, 
the  expected  value  and  the  variance  of  the  contact  number. 
•  r.ixe-  -respective! 


-  r.v  .        j  -  :.: Cp- 

-  tt*S     ":.),     ;  ■   (p-J  (3) 

:-.er  hand,  we  have   a  cne-pc 
segregated  state,  etveer.  different  particles 

boundaries  are  n-  .«., 

Pr{Y  -     Chi  sample  particle  is  oi   component     -  1. 

j  -  i.2 (p-1) 

?r  V.  •  0   the  sample  particle  is  of  component     ■  0, 

(p-1) 

means  t      e  expected  value  and  the  variance  act 

ited  state  are  both  zero,  i.e.. 


EfY.]  =0,  ;  =  1,2 ?-l)         (6) 

s   3 

[YJ  =  0,  :  =  :.: -1)        (7) 

s   j 

Fron  eqns.  (2)  and  (6),  ve  find  that  : 

on  component  .-..    changes  fron  zero  in  the  zc~             I          -c  statr 

n*X.  in  the  complete!  ec  state.   Akac  et  al.  [2]    nave  usee  :  i 

to  define  the  degree  cf  mixedness,  K,  for  a      eneous  binary  sys: 


Cl(0)  "  £s[CI(0)3 

■•  '  -  rc   i  -  ~  re   l  c  * 

where  C.  ,~.    is  the  mean  contact  number  of  component  .-. .  estimated  from 
sarr.oles  of  the  mixture. 


C.   Estimation  of  the  Mean  Contact  Number  by  Spot  Sampling 
.en  a  spot  sample  of  site  n  is  drawn  from  Che 

the  concentration  of  particles  of  component  A   is 

x,  =  -^  ,  j  -  0,1,2 ;?-:       (9) 

where  a.  denotes  the  number  cf  particles  of  componer.' 
sample.   The  mean  contact  number  is,         - , 

cj(0j  -  n*Xj,  ;  -  :,: ■-.  (10) 

Suppose  that  k  spot  samples  of  an  identical      n  are  obtained 
from  the  completely  random  mixture.   The  contact  num.:-:  i  by 

.  ;les  of  component  A.  per  key  par 

[CJ(0))l  -  [n*x    .  L  -  1,2, 
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Equations      -      and      L5      Imply   that    the   estimator,    C.,0»,    is  bi  it 


is   consistent,    since 
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E.   Precision  of  che  Population  Mean  Cor.      .umber  «. 

The  precision  of  the  escitr-ator ,  C'  0  .     ei  population 

mean  c  .s  decenninr     ?ugh  evalua:  .  the  variance 

tribution,  uh;  (see  Appendix  B) 
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3.   DEGREE  DF  MIXEDNESS     il  ON  THE  SAMPLE  VARIANCE  AN!  ESTIMAT1    )F 
ECISION  DF  THE  DEGREE  DF  Ml     ~-SS 


Over  30  different  criteria  for  the  :=::ee  of  mix  have  :een 

reviewed  anc  s  jrr-.arizec  by  Far.  e:  al.  [1].   Here  we  are  ccr.cer:.  . 

wit]  the  cs:ir.a:or  of  the  population  variance  of  an  ind  ent 

from  soot  samples.  ~7,  which  has  been  frequent]  j       . 

J  .... 

degree  of  mixedness.   Suppose  that  a.  denotes  the  number  of  par:.  L< 

conponent  A.  In  a  spot  sample  of  size  n.   Since  a  spot  sar.rle  is 

be  coi         Lxed,  "-.is  random  variable  a  '?se.         -- 

distr.    .  a    i".  the  completely  mixed  state  if  ve  cr.oose  A_.  as  at 

and  lump  all  the  remaining   components  as  the  other  (pseudc 

Thus,  we  have 

E  ra.]  -nX,,         :  =  l,2,...,(p-l  23) 
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=  nX^         .         j  -  1,2 Cp-1 
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listril       can  be  approx.  _.-;:.- 

butioi            •  [7].   Therefor  .ng  ap;      are  tq    .  :-.s 

can  be  obta.              roper ty  of  the  distl  I  . -ince  c: 
normal  distrJ  itJ r.    (see  Appendix  C)  . 
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0   1       p 

Estimates  o:  Lei  di: 

can  be  0'  '  follm 

ng   relationships  are  obtaine  g  each  senple 

mean,  y  ,  to  its  i        do,  eqn.  (16): 
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p 

J  1 


J   -  -.P  (20) 


•  re 


P 


i-0   a. 


i   w«   have 


-5     - 
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- 


.IS 
J  J 


- 

«2.J=2 


,    — 

(1+1 


ins.  2),     -e   ob: 


- er   of    spot  .    -  • 


-    10 


j    ■   1 , 2 , . . . ,  p 
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I    e    ?ara=e:er,     "    ,    car.   be   obtained    from   eqns.     (21]    ar.d     '.2-      as 
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i-l  i 
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-.   PRECIS]    )F  THE  ESTIMATION  Or  MEAN  :::;TACT  NUMBER 

ve      :on8ider  the  precision  cf  the  estimator,  C  .fl>,  eqr..  [6  . 
-  mixture  Ln  an  incompletely  mixed  state.   This  rretis- 
from  the  Lstribution  giver,  by  eqr.  (12). 

The  concentration  of  A.  particles,  >:  ,  in  a  spot  sa: 

j 

expressed  as  the  ratio  of    the  numbei     ..  par:.  Les,  -  .  CO  thl 
sample  sice,  n,  i.e., 


1    ■  - 


eqns.  (16)  and  (17)  spect. 


■i"j'-i«i 


P 


-  11 


ap.«: 


J 


. 

h 


•  •qr.s.  CB)i  (27)  and  (28)   Che  cxpe  »4/ft\]i  and  the 


j(0) 


ill  be  derived  as    [see  Appen-     : 


- 


o. 


j(0)J 


and 


VlCj(C 


• 


*     :.-: 


-- 


(30) 


and    (30),    ve  e    Stand*]  -S.Z 


■  ■ 


fC    ._     . .  -ace   as 


0)' 


-  12 


"-  C 


:;(0) 


: 


-   .   -     -       -     -.7 


:   i-]  •-_--: 


1 


_ 


+  (n-2)x(i+r)(e  +eQ)  -  n+r+i  -  l+i  . 

By  setting  ?  =  1,  this  reduces  r.o  a  known  expressicr.  for  -^  binax       re, 
.,  [3] 

r    (t*+ 

p  c  r  •'  r '     =  . ~  ■  -  _  •  -  - 

1(0)-   n(n-l  ■:  r+2   --;-.-.   

I* 

2n  -       .-,  -  -  :-    -  -  -  (32 

Equation  .        samples,  Che      _ve 

standard  error  to  be  within  a  specif.        .  I.e., 

k  -    ( s [2        -      - 

n(n-l)(x+2 

o  : 

7 

2nr+3)en9      -     n-2)i    l+i  -         *  (n+H  H  (33) 

•    '  rrespor. ..         to    tht   con  .  s.     (30)    t 

;ns.  (10)  and  (11  ,  r*tp«<  ( 

(33)  becc; 

i  ♦  (  Lx. 

,      0     •  0 

■ ' j-1 

ft<  R.S.E.  (Cj(()))Z 

Accor  .  eqn.  .  •  .  .   1  .  ressed  as 
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5.   E5         OF  SAMPL  PRECIS: 

Various  de:      ns   of  the  degree       edness  based  on  the 

sample  variance  have  be        -ed  and  s  in  et  al. 

Here  we  are  concerned  u  .e  estimator  cf  the  population  variance 

of  an  individual  component,  r  ,  which  has  been  frequently  employee 

■  degree  of  mixedness.   Suppose        spot  samples  of 

are  take     random  from  a  ~  I  •  incomplete      /.ed 

•  .  -  C  the  rancc  are  d is- 

I 

buted  according      e  Dirichlet-multinocial  c.        on.   Sample 
mea  .  .   ecr.s.  (16)  and  (17)  ,  are 

repeated  here  for  convenience. 


.1  -  nX.,  -  1,2 (16) 


i]"°] 


L-fj).     J  -  :.: ("> 


:ore,  the  expected  value  and  variance  of  sample  variance  C     poiMBt 
are  [6] 


E[S:;  -  (36) 
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and 


v[s2]  -  f-_   -  ^   ' 

J  -  --   I 


for  it  >  1 


where  -..  denotes  the  4th  c^::^l  moment       about  E[Y  ],  [see  A        D], 

J 


-  E  (Y  -  E[Y  ])* 

~ "  j  J 


Lat  Ive   Bta     lard  error   ol       .  .    R.S.E.  .is 
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Y   and  Y   h..  lomial  J.ist:           •>  parame: 
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(3)   A  set  of  multinomial  random  variables  ,  each  vitfa       -ters 

n,  Z   Z        Z        is  then  generated  [12,13]. 

Steps  (2)  and  (3)  are  repeated  until  k  realizations  of  t 

random  variables  (Y  ,  Y0)  have  been  selected. 
(5)   The  concentration  of  each  component  is  calculated  ace. 

eqn.  (26),  for  each  spot  sample.   The  mean  contact  (,ft. 


■  / 
'2(0) 


and  C   n- ,  are  then  estimated  by  means  of  eqn. 


(6)   Steps  (2;  to  (5)  are  repeats     tines,  and  the  Bean,  ?'  nv, 

and  the  variance,  V[C!  ,..],  are  computer  :  rs , 
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Substituting  eqns.    (D-6      through   (D-9)    into   eqn.      D-1C    ,   we   final] 
the  ceiailec   expression   sf    the   fourth   central   moment   of  tl  ilet- 

nultinomial    cis:r: tioti    for    component   .-. .  . 
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Statistical  analysis  r.as  beer,  a  major  tool  of  solids  nixing  investiga- 
tions because  of  the  stochastic  nature  of  nixing  processes.   Over  thirty 
different  criteria  for  the  degree  of  r.ixedness  have  beer,  reviewed  and 
E  .-r.arized  by  Far.  et  al.  [2).   Tr.e  relationships  among  so~e  cf    the  criteria 
have  been  studied  and  tabulated  by  Far.  and  Wang  [3].   Table  1  lists  these 
criteria  and  investigators  who  have  proposed  then.   Most  cf  the  criteria 
are  based  on  the  composition  analysis  of  equi— sized  particle  samples  of  a 
nary  mixture.   However,  sclid  particles  tc  be  blended  often  consist  not 
only  of  particles  cf  varying  size  and  density  but  also  of  widely  d:::rr- 
chemical  ar.:  physical  properties.   Tr.us,  a  study  cf  mixing  indices  for 
D -1  ti component  heterogeneous  solid  particles  is  of  practical  sign:..     e. 
A  convenient  scheme  for  clas.-    tag  particle  -     ^.  should  be  es 

facilitate  the  investigation  cf  sclids  c.        Usually,  par: 
systems  are  classified  as  homogeneous  or  heterogeneous  vor  dispersed) 
systems.   Particles  of  a  homogeneous  system  have  the  sane-  physical  properties 
(e.g.  den;.-  ,  size,  shape,  etc.)  and  are  distinguisl  the 

:  ics  cf  tracers,  such  as  color  and  radioac: .      ~:.ich  dc  r. 
course  of  the  mixing  process.   He:  « 

physic.  _  .:h  do  not  interact  nil   the  opt:     .  of  a 

are  considered  as  tracers.   On  the  other  hand,  particl  .eous 

system  have  different  physical  properties. 

The  extent  of        .at  ion  in  the  c:        .sties  of  a  mixture  may 
be  •  ::.s  of  v..  Lght,  volume  or  number  fraction  of  each  con:'.  ::- 

:e  [4,5).   The  three  bases  yield   an  identical  r>      I  or 
a  homogeneous  system.   However,  sampl .      counting  the  number  of  particles 
is         t  and  impractical  for  a  fine  particle  or  powder  system.   Measure* 
lume  fraction  is  not  easy,  because  the  bulk  volume  of  sc 
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particles  consists  of  the  true  volume  of  partnles  and  t  ze, 

i.e.,  the  bulk  vclume  of  each  cor.s:::.er. :  is  net  at-      .   Since  the  practi- 
cal  difficulties  of  measuring  the  •       ix    a  nixer  exist,  the  zear.  ar.t 
variar.ee  of  the  weight  fratticr.  are  chosen  here  as  unit      easure 

I  te  literature  or.  the  assessaer.:  : :"  mixture  quality  is    g  ry 

scarce.   This  is  especially  true  for  mul  ticomponer.t  i      es.   Gayle  et  al. 
[6]  used  chi-scuare  as  a  criterion  instead  of  conventional  defil 

.  tng  index,  which  utilize  the  sample  variance  ior  sta:.:„r^  deviation 
Granules  used  in  their  experiments  were  sir.ilar  m  surface  cl  I         --S, 
sr.ape  ant  density  but  different  in  eclor  only.   Luslik  [7]  propose^ 
numerical  homogeneity  index  to  express  q..  y  varying  degrees 

homogeneity,  nar.o     the  negative  logarithr.  cf    the  sample  w< 
to  obtain  a  standard  deviation  o:    1*. .   Heri.ey  et       .■   jse       con- 
cept to  ca  ty  of  ad      ingre_. 

hey  also  assumed  tha:  |    components  z       ■  to  the 

same  particle  size  level  prior  t    I  :ng. 

Ifl  the  present  at    .  some  approaches  to  assess  the  c  a  i:>.- 

■VC  been  propos  roxim^  cneous 

derived. 

2.   51   .  X SIS  OF  BINARY  MZXTU1 

Thi     tsaioni  of  some  statistical  quantities  associate.       a 

are  s      zed  below  [2): 
A  sample  variance  is  defined  as 

'  -  «  Z   (xli  '  (1> 

-1    1J 


An  ur.biaso-  sa-        ance  is  defined  as 


2 


ere 


■  '  s . 

x.   ■       or.  of  compc:  th  saa: 

-  sarrie  Dean  which  is  -   I   >. 

e  cal  :ncr.:       -qual  tc  that  cal: 

:re.   The  popular:         or 

k2]  -  (3) 

e   see:  rd    deviation   or   variance 

eiy    rand-  ed    state   and    I 

segregated    state   are   cfter.   used   as    the  boundary   values 

aoong   a    group   of   samj . 
drav  the   co-pletcly    rancor    state    is    ce: 

[10)   as 


-:  -2 

- 

-...-.'_■  (4) 


- 


V  V 

I 


re 


.    «   ve::     :  and    2,    respective. 

-         (Px  +  P2  -    1), 

-     .       El    of    compone- •  and    2,    rttp€£td    >.     . 

i        - 

s, .     -  standard  d      Lona  cf  7   •  of  components  1 

anc      ■ spective 

•  -     C  of  I  sar.rie. 

.e  is  ..  (4)  is  reduced  to 


-  5 
P,*2    - 

If  each  component  consists  of  particles  c:  Che  s.v      :.;le  weight)  w, 
(w,  =  w_  =  w) ,  ecr..  (5)  car.  be  written  as 

2     V  '   T, 

i  w  :  : 
=  n  p:p2  (6) 

h  is  the  original  formula  of  Lacey  [11]  where  N  is  I       ei  c:  par:  .- 
ties  in  a  spct  sample. 

~er  the  assumptic:     t  the  pro:  el  Llity  of  sar 
the  given  component  is  proportional  tc  its  volume  fraction,  t.ne  variance 
for  a  mixture  in  the  completely  segregated  state  is  derived  here  as  .'see 
Appen- . . 

o     ■  —  CI— I  I  P  -p  "(1  -  —  P  ) 

-  -■-   !-:•  .  :  :■  - .  "  1  -  —  P2)  (7) 

vheri     end  .  .  denote  trie  densities  of  components  1  and  .  .  .v, 

and  z    thl  the  components  have  the  same  particle  de- 

(p,  ■  p,  ■  e)i 

-  P  P  (8) 

can  be  seen  that  th-     .ance  of  a  mixture  in  the  completely  segre- 
gated state      (dependent  of  sair.;  .ded  that  sampling 
taken  at  the  interface  between  the  components  and  tha:      variance 

■  concentration  (.         ion)  :  d  sta: 

la  inversely  proportional  to  sample  sise.   Thl  : en  most  free       used 
.ces  for  binar      ires  are  the*     Table 


V   -   t 


A 


a  r  r «. 


geneo 

■ 

tin' 

0  (9) 

P 

iough  spc 

poslton,    such  a   :  Lng 

op- 

state  thai  is  d it  ideal  homogeneity.  It  i-  ed  by 

of  sf.   '     ..  ■  e  components  from 

.lure  ,  even  though  the  c«.  I  be 

tdenti  e.  For  a  homogeneous  binary  system,  particles  in  a  t| 

sample  of  ■       are  distributed  according  tO  I  ial  difttt  .    •   . 

Pr|V  V  ^Tri"la-r/2  (10> 

where   V      i<-  ble    that,    denotes    the   DUiabei  .cs   of   cor 

ponent    lira.'  on, 

th<  .  s 

o      «   Var    [x, ] 

r  1 

-   Var    [J  ) 

which    is    eqr..     (6)  . 

In    pr  Lng   pro..  • 

C    because    the   pro  or   b'  :<y 

d  i : 


thai    the   probability    distributions    of    the   components   curing  .    are 

Down.      The   probability     :f    sampling    a   particle   of   a    given    component    is 

proportional    cc    the  volume    fracrior.   cf    rhar    corponer.:        Dnder    che 

•    -jsis   cf   mixtures    In    Che   completely   randor.  scare,    parricles    ir.  a 
sample   car.   be  assumed    to   be   apprcximarelv    ciFtri:  according    re    rhe 

binomial    distribution  with    parameters    N      s.r\i    q^  ,    where   K      is    the    •  .  ed 

b  er  of  pamrles  ir.  a  sample  cf  constant  weight,  W,  and  c.  is  the  vcl 
fraction  of  conponenr  1.  Tne  variance  cf  rhe  weight  fraction  cf  coupons 
1    or    2    is    derived   as    [c  :...:■    B] 

pl°2         12  f2  1 
Since,    ir.    practice,    solid   particles    to   be   r.ixec   are   usual]  isured 

sar_  sns ions,    this  distri  seems    I  a  good  ■ 

■  -res.       Notice    that    ear..       11       reduces    to    eqr. .     ('  hocogen- 

eous    system. 

3.  [CES  FOR  H  IICOMPONE     ::TUKES 

Th     pose  c:  •     section  is  to  develop  mixing  li 
cor  ture.   We  SSSumt  throughout  this  section 

:  a  mul 
(by  weight  .corded.  •  sampl  OS, 

an  indication  of  the  degree  of  nixednei  •  v.  - 1 1  be  calculated.       ods  of 
sampling,  determination  of  sample  lise,  es  dr.. 

on  the  degree  of  pri  I  cost  and  practical, 

of  sampling,  and  possit         al  constraints  on  the  sampling  process. 


"Pseuc 
Suppose  that  a       -poner.     curt  o  .... 

aogeneity  of  the  key  component,  say,  component    | ne 
e  is  con        co  com  1  and  another  (psuedo-j  com- 

reaair..:..   . -.    •  onents .  is 

pre:  ..  prop       of  (1-P-).   The  p.-.  •  and 

de:  are  give  as 

-        -:-v 

—  -  --  -  ...  - 

v. 

- 

U-Pj) 

m      

o     p 

L2  '  k 

.    cospc  :  e  treated  as  conponent  2  of  I  bit         re  as 

de^       r.  sec:       ?:r  Instance,  I  Lndustry, 

us  u  ent  and  several  diluents  (lactose,  stare  . 

.cose,  etc.)  are  ir.  a  dosage;  ther<:       the  homogeneity  ' 

key  conpon<      ■  the  ca  ir.  concern.   The  r.  I  .ie 

be  dtl     .e:  as  described  ir.  the  7receding  section  UJ, 
teh  based  on  the  concept  G  -re 

0  be  a;       tc  each  component  in  a  oulticotr.poner  '         to  define  an 

each  co~  •:  thod  is  acre  useful 

.  cate  that  a 
for  II    .oaponer.t,  i.e.,  sone 
ts  nay   be  requ:-         mil      a  nuch  greater  e>  other 

(1). 


-  c 

3.2   ?c:led  '.ariance 

Zr.e    saccle  variance  cf  the  concentration  ir.  a  -  iti       ent  mixture 

is  defined  as  follows: 

k    ■ 
9    --^r-  Z   '  Z   (>:.  .  -  x.  )V  (14 

m— 1  _  l     ij    i .   i 

where 

k  =  number  of  cor.pcnents, 

c  -  Dumber  of  spct  sanples, 

>:  =  fraction  cf  component  i  In  tbe  j— th  spct  sample 

>:  .  =  sar.ple  near,  of  component  i 

=  -  :  x±i 

When  k  ■  2,  eqn.   14   car.  be  reduced  re  eqn.  (2  •  -      expresses  :      ::- 
ance  of  i  ture.   S.     ly,  we  may  express  the-  variance!  for 

uulticor.ponent  mixture!  in  the  completely  mixed  and  segre,  . 
respe         as 

c:  -  :  k.1  P]  (1  -  —  V  (15) 


and 


1-1  111 


.    -  ompon«.  res 

as   sir.  to   those   for  binar  Table   1.      Pol     the  homogenous 

Par'  tta,    eqr.s.    (15)    and    (16)    can  be   red.         .    reaped  .    :o 

(1-PJ 
1 


0    :-. 


- 


.  -e  have 


,  - 


Ac  ne  relat.     .  ps  bet-     :hese  I        -ices  sh. - 

Table  1  ere  iden:ical  to  those  for  jres  [3].  -   - 

plot  the  value?  of  the  nim  s  aga: 

.ind  200 

car.  be  a<  related  I 

6 

a??-                be  more  depender.  ers.   However, 

the  variat.       Eht  e>:t.        xedne^  :ne 

regated  state  tc  the  c  .xec  state. 
.e  Covariance  M 

inalysil  is  an  effectivt  tool  f.  :  icompone 

:  observing  the  c.  ion  of  a 

of  a  k-cor.por-  erned 

I  es  c:  (k-1)  ve        .ec  rancor  respcr.se  v  |  e:-  .   Nrtice 

oncer         need  tc  be  consider^  .en  sample? 

are  o":            a  (k-l.'-varia  te  popula*        \eac  of  ^  riate  popu- 
lation, tht  data  can  be  represented  as 


.co- 
sample  T\v 


(k-1) 


(k-l)l 


(k- 


V  -  10a 


. 


"■.DFX  K 


B 


.  1   Comparison  of  mixing  indices  based  on  the  2nd  approach 
(pooled  variance)  for  sample  6ize  of  20. 
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Comparison  of  mixing  indices  based  on  the  2  nd  approach 
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le   near,   cf    the    l-th    component    is    defined    as 


:•:         =   —      -       x  .  . 

■     =       ij 

The  sins  of  squares  and  cross  products  are  defined,  respectively,  as 

;  ,     -  >* 

=      _       ;  x  .  .    -    x  .     ) 
j-1         «  l< 


2        1 
j-1         "  j=l         J 


• 


and 


"v."       ■      ■  '(^-: 


r. 


r. 


-    x . .  x      -  -  (  :    x . . )  (  :    i 
J-3  c     j-1     «     j-1 


The  estitaated   variances    of   x     and   co .  oe   of   x     and   x     are,    r< 


[x.]    -   si 


SS 


m-1 


1-1,    -.    ....    [k-3 


-■ 


and 


cov. .-_.     ;  -  su 


s? 


All 


.  i 

.    :he  sample   \  jriance  matrix    is 


.  (V-- 


- 


s  - 


.. 


•-: 


Sl(k-1)     .  ,-.  S(k-DU-D 

a  com:  re  are 

expec:         randc                    spot  ;>ns 

the      re.   Let  us  de:       random  vec:  


r.er. r.s,  vnere 


'J 


.  C  ■_    r. 
j 


• 


par- 


:ne  spot 


randc      .ables. 

1    - 


\ ,  "  \  ] 


SU p 

n.  ln,!...iL  I  rl  r2 


! 


i-i  a 


- 


P   -     -...-7,-1 

...... 

«  .  -  .  ture  (for  1 

homogeneous  system) 


-   I 
Note    that    the  ::_!  ttnc~ial    distribution    is    a    ger.eraliza: .  _      -  .al 

distribution.       According    tC    this    pr       l    ty    density    function,    ecr..       25), 

the  ir.eans  ,    variances    and    ccvariances    of    the   rai  .es,    V.,    Y«,     ■•-,    V. 


,re   given,    respectively,    by:       [14,15] 


i  .   ;  ■  w> 


VarlY. ]   =  NP,(1-F.) , 

Cov[Y..  ,    T   ]    =   -HP.P. 

l        i.  _    - 


21 
C27 

-- 


or  equivalent!;  . 


Elx,]    ■    Ehf] 


=    V 


I   =  Var[- 

Is 

P.(l-P.) 

l  1 


(30) 


Gov 


■ 
,   xt]   =  Cov[-jj-   .  ■ 


P.Pr 

i    x 


(31) 


Therefore,    t  ..ce-covariance   ma:  the    C< 

can  be    expressed   as 


PjUpPj) 

-P1P2 


'P      P 

-1) 


-P     P 
P2(l-P2) 


-  I       P 

?(k-l) 


- 


:'(k- 


I 


P(k-l)(1  1)} 


' 


segrege^      :ate,  each  different  cospc: 
separated  a  |  pe.   I         ate,  we  have  | 


end 


Var[xJ  - 


■  -PiPi 


(33) 


:ance  cat  • icocpone: 

can  be  expressed  as 
.   .   . 
P2(l-        ... 


l~Pl> 


?, 


i 


p 
:*(k-i) 


_  i 


(l-p    ) 


(35) 


e  variar.ce-cevariance  r      .  £,  characterises  the  degree  of 
the  dispers  coraponer.:.   Xhue(  Chf 

.        the  extern 

riencf  ir.  the  u:  .alysis.   Therefore, 

tnd  0-  by  S I  ,   I  I  ,  and      ,  respectively.  -   car.  also 

Dent  B  -is,  we  can  extt 

ex 

2 

-  - 

(36) 


0-0 

m  L 

:   : 

°o  * 


:ey  [16]  for  a  I ure  to  a  multicor 

ng 


- 


-s 


are  t  the  expressions 

are  reduced  tc  those  :  "ures. 
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According  t:  eqns.  ( 22y  and  (25;,  we  have 
cr.: 


The  relationships  among  the  mixing       es  tabulated  in  Tatle  1  are  evalua- 

ir.  Table  2.   Figures  3  and  6  plot  the  values  oi  I   nine  n        -ices 

against  that  of  M-  for  the  sar.rles  of  sizes  20  and  20C  ,  respectively.   A 

mixil   Lndex  should  be  sensitive  to  changes  ir.  tl  c:  a 

mixture  ir.  order  to  be  useful  ir.  following  the  progress  of  mixing. 

Vaier.-     L7]  pointed  out  that  Lacey'f  index,  ecn.  (3"),  was  very 

sen:         changes  near  the  completely  segregated  state  and  rela: 

insensitive  during  the  process  of  nixing  near  the  ccr.pletel;         :  :e. 

Their  more  complicated  index,  !'._,  was  thus  devised  tc  k  :er 

sensitivity . 

:.ed  by  the  determinant  :e- 

covariance  matrix  for  mul  t  iconpone     ..tures  can  be  gent: 

.  where  I      and  Z      are  the  var iance-covariar.ee  mi       c:  w 
— r     — s 

proportion  in  t.      pletely  ^r.ixed  and  completely  segregated  stau  . 
respectively.  The  development  c:   _r„  an^  -     Frc:St 

only  with  homogeneous  systems.  Co:  eterogeneous  case.  ut 

probability  of  selecting  a  particle  of  a  given  compont      proportional  to 

volur.<  :hat  conponent  and  the  probabil.  I  ng  a  particle 

ne  compone:.-        tical  everyvhi      a  coaplet<     andoc      re. 
Under  these  assumptions,  particles  take:      a  spot  sample  are  die:     ed 
approximately  according  to  the  multinomial  distribute         e  parameters, 

where  N   is  the  expected  number  of  particles  in  a 
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Fig.  3  Comparison  of  r.ixinp  indices  based  on  the  3rd  approach 
(  the  detcrrirj.int:  of  the  covariance  matrix)  for  sample 
size  of  20. 
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.    4      Comparison   of   mixing    indices   based    on   the    3rd   approach 
(    the   determinant   of   the   covariance   mat:  r   san-. 

size   of   200. 


:r»e 
componc: 

.mces  and   c:  -  •  .ons   cf    components    in  a 

joporu  e,    ha  cen;  >re 

re  Appenc | 

•■      :  -     -  p<>.  0) 

and 

- 
Ccv  .     :  ,     |  :         . 

,2 '  (41) 

.  ishes  .nee   m  .re 

•       c,    I    .  cht   variances   and  covariances 

■     I    |    components    in   c 

the  segregstec  are   givi     .    respectively,    by    [see  Appendix  D) 

]  -  £■  (l-l     *  l     -  •  -  -1-  p.   .     i  -  L,   : -: 

:       .    .      l]-Pj         L-         -         . 

.  I    -  1,   : k  and   1  t  I  (*3) 

es    the   var iance-covar iance   matrix    for   a  .ompor.  :re 

te,    £    •      Not       e    that    the   variances   and   co- 
•rogeneous    I  I  reduced    to    those    for   a  homogeneous 

px  .  o2-   .    .    .   .pk 

■  ■  - 

... 


.     17 


1 .  -  .      Sa-ple  Mean  of    ar.   Individual    "omponent 

Ar'-ac    et   al.     '!£]    have   proposed    a   mixing   index  base,:  the    zz?.- 

b  er  as 

C  -  C 

s 


c    -  c 

r  s 


r 
where  C  is  ar.  estimator  of  the  population  contact  r.  a  -;>::ure 

ar.  .-  .irary  state,  C^  arid  C   are  near,  contact  nurDe::  of  th<  Bixturei  I 
Lately  r.ixed  state  and  completely  segregated  Bt     respect  J   ly. 
index.  :..  .  ,  eqn.   --  ,  taK.es  the  value  of  C  In  the  completely 
segregated  state,  a  value  betveen  0  and  1  li  ...  in  .xed  state, 

an.        e  COnplei         d  state.   Recently,  Shindo  et  al.  :st: 

Lnomial  mod<_      ^escribe  zhe   mixture  in  ar.  incomplete         state 
Che  parameter,  t,  vhich  is  unique,  car,  be  employed  to 
Lng  index 

".:  "  —  (45) 

shown  that  t'r. .       Ltion  can  be  extender  | 
Bixing  index  as  [20] 

Kll  "  771  CW 

The  relative  standard  error  of  the  mean  contact  number  estimator  has  been 
sh,  -      be  sn.        ,n  that  c:  |    variance  estimator  of  spot  samples  120] 

be  defined  by  using  the  sample  mean 
each  ind.     I  compon.        tJd  of  -!,ich  is  con. 

used  to  define  the  i      inaex.   For  instance,  I  g   expression  U 

be  used  as  a  r      index  for  a  bina-      :ure; 


LI 


. — i . (47 

The   de:  equa:  .1    sur  e   abs  ences 

ifl   sar  s   and    th<  conct  .on  among 

the-  >-s    froc  a  cota:  segregated   state    (See 

hi   populi.  .       •    can  be 

le  near..      The   co:  re 

DUghoul  %  equc.  I  ,    and 

therefore,    the  :he   pe:  I    accc: 

•    e   co-.pletely    segregated    state, 

ecr. .      47),         -  takta    ;--  between   C   and    1 

also  be  ex: 
.t    system  as 


:  f    :       :•:      -         ^  p. 
■  : ,   ■"  ~x- - 
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±  Ci-p  )  +  (1  -■£■  p  )] 

.     i 


4. 

Tc  was    choser. 

to   c  use    (  o posed  ices    for 

es.      Iht   values   of    the    degree   ol  ss    ev.  d   according    to 
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Appendix  5.  Equation  (11) 
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Append  f  the  Absolute  Differences  between  the  Staple  Concentrations 
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CHAPTER  VI 
CONCLUSIONS  AND  RECOMMENDATIONS 


This  chapter  summarizes  significant  conclusions  of  this  work  and 
gives  several  recommendations  for  further  studies. 
1.  CONCLUSIONS 

(1)  The  concept  of  the  contact  number  to  estimate  the  mixing  index 
is  successfully  extended  to  multicomponent  mixtures.   Such  an 
application  has  been  employed  only  for  the  mixing  of  binary 
mixtures  111.   The  multinomial  model  is  proposed  to  describe  a 
multicomponent  mixture  in  the  completely  mixed  state.   H      ~al 
experiments  have  been  performed  to  verify  the  theoretical      i s- 
ions  of  the  distribution  of  the  contact  number.   The  mean  contact 
numbers  and  their  precision  from  the  computer  simulations  of 
mixtures  with  regular  packing  arrangement  are  reported. 

(2)  A  theoretical  model,  namely,  the  Dirichlet-multinomial  model,  - 
derived  to  describe  a  multicomponent  mixture  in  an  incompl 

mixed  state.   This  model  contains  parameter   that  defines  uniquely 

based  on  the  contact  nu:  ■  ...is 

parameter,  I  ,  is  derived  and  the  precision  of  the  estimated  po;  - 
lation  mean  contact  number  is  assessed.   Computer  simulations  a: 
included  to  demonstrate  the  effi  he  derived  modi 

I  andard  error  o:     MAO  contact  nu 
has  been  shown  rstima. 

of  spot  sample*.  the  precision  index 

based  on  the  mean  contact  number  estimator  is  su.  at 

CO  based  on  the  variance  estimate 


- 

(3)  Literature      ■  assessment  of  mixture  cu.  scarce,  as 

Ls  topic  it     acks  generally  accepted  principles,  in  spite  of 
practical  significance.   This  ls  especia      ue  for 
component  solids  mixing.   Several  approaches  to  assess  the  homo- 
geneity of  a  multicomponent  solids  r.ixture  have  been  proposed 
Chapter  V.   An  approximate  distribution  of  a  hetero^ 
in  the  corcplet'  also  reported. 

proposed  here  are  estimated  on  the  I      of  the 

(a)  pseudo  binar  re  conct 

(b)  pooled   variance  of    the  whol' 

(c)  determinant   of    the   sample   covariance  matrix,    or 

(d)  sample  mean  of   each   component    in    the  mixture. 

(1)      The  non- ideal  mixing  behavior   of   cohesionless   solid   particles    in   a 
•  izontal   drur  r>een  mode  a   continue 

process    [2].      The  Markov   chain  may   also  be  applied    to   Simula- 
the  non-ideal  mixing  behavior   of   solids   particles    in  a   dr 
Lai  and    Fan    [3]    proposed   a   discrete   steady-stat  kov   chain 

model    for    the  mixing   process   of   a  multicomponent   homogenec 
particle   sy  D  a  motionless  mixer.      Wang   and   Fan    [«]    also 

proposed  a  Markov   chain  model    for   axial    I  :tion  of   bir. 

solid   particles.      It    is   recommended    that    this  model  ma-  •  s- 

••.component   heterogeneous   particles  blended   by   a 
motionless  mixer. 

(3)      Other  »f    stochastic  models,    c.  ne  birth-death  model 

(5,6]   and    the   restricted    random  walk  model    [7,8".  be   employed 

to   simulate   solids  ■  and   demixir. 


(4)  Multivariate  analysis  is  an  effective  tool  I  ilti- 

component  solids  nixing  problems  [9,10,11].   Since  tl 
are  statistically  dependent  among  themselves,  v  jle 

out  one  from  the  others  for  treal     In  bulk  saapling        Lti- 
component  heterogeneous  i:.i>::ure.   The  vai  be  cc 

jointly.   Thus,  the  following  multivariate  statistical  methods 
may  be  exploited  to  analyze  sarpling  results: 

(a)  multivariate  analyses  of  variance  a: 

(b)  multiple,  partial,  ar.d  canonical  cor  hniqui 

(c)  discriminant  and  classification  analys.  . 

(d)  factor  analysis  and  canonical  ana 

(5)  Nonparametric  statistical  tests  can  be  :  data  v. 
different  scales  of  measurement,  such  as  nc     . 

or  ratio,  without  knowing  the  distribution  cf  |      pulation  [12), 
The  distributions  of  the  components  during  mixing  are  usua. 
unknown,  and  therefore,  nonparametric  stat;  :iould 

be  effective  tools  to  anal        sampling  i  olids 

nixing.   Applicability  of  these  |     on  solids  mixing  i 
mi        :s  demonstrated  [13].   Natural  .in  be 

omponent  mixtures  and  r. 
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ABSTRACT 

The  concept  of  the  contact  number  is  extended  tc  e  oultic<      at  ioli 

mixture  in  the  completely  mixed  state.   Chapter  3  presents  an  es       cm  of 
the  mean  contact  number  by  measuring  the  concentration  of  particles  by  s: 
sampling.   This  chapter  also  reports  the  mean  contact  numbers  and  their 
precision  determined  from  the  computer  simulation  of  mixtures  with  regular 
packing  arrangements  at  different  concentrations  of  the  key  component.   The 
number  of  spot  samples  necessary  for  a  specified  value  of  the  relative 
standard  error  is  also  derived.   All  equations  derived  here  are  reduced  to 
the  corresponding  expressions  for  two-component  mixtures. 

A  theoretical  model,  namely,  the  Dirichlet-multinomial  mod-  '.  .       i  .  vec 

for  the  distribution  of  (Y.  ,  s.'  ,  7  N  In  a  spot  sample  oi   size    r.  of  a 

12         p 

multicomponent  mixture  ir.  an  incompletely  mixed  state,  where  Y.  denotes  the 
number  of  particles  of  component  A..   This  model  contains  the  parameter .  7". 
that  defines  uniquely-  the  mixing  index  based  on  the  contact  qui        -timate 
of  this  parameter ,  7",   is  derived  and  the  precision  ol  Na- 

tion mean  contact  number  is  assessed.   Computer  simulation  is  included  to 
support  the  effectiveness  of  the  derived  model. 

.terature  on  the  assessment  of  the  quality  of  solids  r.ixtu: 
very  scarce.   Tins  is  especially  true  for  multicomponent  mixtures.   Four 
approaches  to  assess  the  homogen.      :"  a  mu         ent  solids  I      I  have 
been  |      »-d  in  this  study.   The  n      indices       ed  are  estimated  on 
the  basis  of  the 

pseudo  blfU 

.ed    variance   of    th-  <•    svstem, 

(c)      de terminal  '  sample    covariance   ma 

sample   mean   of   each   component    in    •  re. 

The   use    of   various   mixing    indices    for  multicotnpon«  es   has    been 

throtlf  ternarv   mixtures. 


